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Abstract 

Currently,  fuel  cell  technology  plays  an  important  role  in  the  development  of  alternative  energy  converters  for  mobile,  portable  and  stationary 
applications.  With  the  help  of  physical  based  models  of  fuel  cell  systems  and  appropriate  test  benches  it  is  possible  to  design  different  applications 
and  investigate  their  stationary  and  dynamic  behaviour.  The  polymer  electrolyte  membrane  (PEM)  fuel  cell  system  model  includes  gas  humidifier, 
air  and  hydrogen  supply,  current  converter  and  a  detailed  stack  model  incorporating  the  physical  characteristics  of  the  different  layers.  In  particular, 
the  use  of  these  models  together  with  hardware  in  the  loop  (HIL)  capable  test  stands  helps  to  decrease  the  costs  and  accelerate  the  development  of 
fuel  cell  systems.  The  interface  program  provides  fast  data  exchange  between  the  test  bench  and  the  physical  model  of  the  fuel  cell  or  any  other 
systems  in  real  time.  So  the  flexibility  and  efficiency  of  the  test  bench  increase  fundamentally,  because  it  is  possible  to  replace  real  components 
with  their  mathematical  models. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  the  recent  years  there  was  an  increasing  interest  in  fuel 
cell  technology.  In  particular,  the  polymer  electrolyte  membrane 
(PEM)  fuel  cell  has  reached  a  high  development  status.  This 
development  was  mostly  advanced  by  the  automotive  industry, 
because  fuel  cells  are  suitable  to  substitute  the  fossil  fuels  and 
also  to  provide  an  environment-friendly  propulsion.  But  there  is 
also  a  growing  market  for  stationary  fuel  cell  applications,  e.g. 
for  cogeneration  of  heat  and  power  (CHP),  and  as  a  substitute 
for  batteries  in  portable  devices,  e.g.  for  laptops. 

Despite  the  progress  made  in  first  realised  fuel  cell  vehicle 
fleets  and  CHP  units  there  is  still  lot  of  work  to  be  done  to  bring 
fuel  cells  to  the  market.  The  reduction  of  still  high  costs  has  the 
biggest  priority  but  the  build  up  of  hydrogen  infrastructure  is 
not  less  important.  Furthermore,  the  durability  and  reliability  of 
conventional  systems  has  to  be  reached. 

By  means  of  mathematical  modelling  the  development  and 
design  of  fuel  cell  systems  can  be  highly  accelerated.  There 
are  many  three-dimensional  steady- state  PEM  fuel  cell  mod¬ 
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els  [1-3]  described  in  several  papers.  The  results  are  the  local 
concentrations,  the  water  amount  of  the  membrane  and  their 
influence  on  the  current  density  [4,5]. 

In  particular  for  automotive  and  portable  applications,  the 
knowledge  of  dynamic  behaviour  is  of  vital  importance.  There¬ 
fore,  the  PEM  fuel  cell  model  which  is  presented  in  this  work  is 
capable  of  predicting  the  dynamic  fuel  cell  behaviour  depending 
on  the  operating  conditions. 

The  use  of  test  benches,  which  can  be  combined  with  the 
mathematical  models  of  the  components  that  are  not  avail¬ 
able,  opens  up  fully  new  possibilities  in  the  development  of 
fuel  cell  systems.  The  test  benches,  which  are  represented  in 
this  work,  dispose  a  hardware  in  the  loop  (HIL)  interface  that 
provides  fast  data  exchange  between  the  test  system  with  the 
available  hardware  to  the  models  of  those  hardware  which  is  not 
available. 

As  a  result,  the  behaviour  of  a  fuel  cell  system  with 
dynamic  loads  is  calculated  with  respect  to  voltage  drops, 
humidification  of  the  membrane,  stack  temperature  and  the 
influence  of  temperature  on  the  reaction.  Afterwards,  the  effi¬ 
ciency  and  the  maximum  performance  of  the  fuel  cell  stack 
and  of  the  whole  system  are  analysed  under  different  condi¬ 
tions.  Finally,  an  application  example  of  the  HIL  simulation  is 
shown. 
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Nomenclature 

^Cell 

cell  area  (m2) 

Wool 

heat  capacity  (Jkg-1  K-1) 

CD 

double  layer  capacitance  (F) 

j 

current  density  (A  cm-2) 

^cell 

cell  mass  (kg) 

^H2,Mem  hydrogen  transport  in  membrane  (mol  s 

Afit2o,Mem  water  transport  in  membrane  (mol  s- 1 ) 

A),gdl 

molar  flow  gas  diffusion  layer  (mol  s  1 ) 

Wsf2,Mem  nitrogen  transport  in  membrane  (mol  s  ) 

Nq2  ,  Mem  oxygen  transport  in  membrane  (mol  s  - 1 ) 

Pi 

partial  pressure  of  component  i  (Pa) 

Pj\\ 

thermal  power  (W) 

Rm 

electrical  resistance  (£2  cm2) 

^Mem 

ionic  resistance  (Q  cm2) 

T 

cell  temperature  (K) 

Th\v 

air  temperature  (K) 

Wool,  in 

temperature  of  cool  medium  (K) 

UAn 

anode  voltage  loss  (V) 

UC a 

cathode  voltage  loss  (V) 

Well 

cell  voltage  (V) 

Uq 

ohmic  voltage  loss  (V) 

Wool 

cool  medium  flow  (m3  s-1) 

yH2,An 

H2  in  the  catalytic  layer  anode  (molm-3) 

FH2,Ca 

H2  in  the  catalytic  layer  cathode  (molm-3) 

FH20,An 

H2O  in  the  catalytic  layer  anode  (molm-3) 

yH20,Ca 

H2O  in  the  catalytic  layer  cathode  (molm-3) 

>7,FF 

component  i  in  the  flowfields  (molm-3) 

yi,G  DL 

component  i  in  diffusion  layers  (molm-3) 

yN2,An 

N2  in  the  catalytic  layer  anode  (molm-3) 

FN2,Ca 

N2  in  the  catalytic  layer  cathode  (molm-3) 

y02,  An 

O2  in  the  catalytic  layer  anode  (molm-3) 

y02,Ca 

O2  in  the  catalytic  layer  cathode  (molm-3) 

O'Ca 

cathode  charge  transfer  coefficient 

^Mem 

swelling  factor  of  the  membrane 

£i 

porosity  of  medium  i 

X 

stoichiometry  of  gas  supply 

2.  Fuel  cell  system  model 

The  dynamic  fuel  cell  system  model  (Fig.  1)  includes 
hydrogen  and  air  supply,  humidifier,  consumer  load,  electrical 
energy  storage,  cooling  system  and  a  detailed  stack  model.  The 
model  has  been  implemented  in  MATLAB/Simulink  software 
and  the  simulation  can  run  with  discrete  or  continuous  time 
steps. 

The  air  and  hydrogen  are  humidified  to  the  set  values  before 
they  reach  the  stack.  The  pressure  and  the  flow  rate  are  con¬ 
trolled  by  the  check  valve  and  the  compressor.  The  cooling 
system  can  be  based  on  a  liquid  or  gaseous  medium.  If  the 
current  or  the  power  for  the  consumer  load  is  set,  the  volt¬ 
ages  of  the  stack  and  of  the  electrical  energy  storage  are  cal¬ 
culated.  With  the  help  of  the  regulation/control  system  it  is 
possible  to  test  a  lot  of  different  operation  strategies  and  the 


power  split-up  between  the  stack  and  the  energy  storage.  A 
further  advantage  is  that  the  system  model  can  be  changed  eas¬ 
ily  because  it  is  built  out  of  different  autonomous  operating 
blocks. 

In  the  following  the  description  of  the  stack  model  is  given; 
the  description  of  the  entire  fuel  cell  system  is  not  the  subject  of 
this  work. 

2.7.  Stack  model 

In  the  one-dimensional  stack  model  the  different  layers 
(Fig.  2)  are  implemented  with  their  different  physical  and 
electrochemical  characteristics.  In  the  diffusion  layer  and  the 
catalyst  layer  the  transport  mechanism  of  hydrogen,  nitro¬ 
gen,  oxygen,  water  vapour  and  liquid  water  is  included.  The 
water  transport  through  the  membrane  incorporates  the  effects 
of  concentration  gradients  and  the  electro-osmotic  drag.  The 
nitrogen  and  hydrogen  transport  through  the  membrane  is  also 
included. 

After  calculation  of  partial  pressures  of  the  gaseous  species, 
the  water  amount  in  the  different  layers,  the  polarisation  losses 
at  anode  and  cathode  and  the  membrane  voltage  loss  are  calcu¬ 
lated.  The  polarisation  losses  are  calculated  under  consideration 
of  the  double  layer  capacitance  at  the  electrodes.  The  outcome 
of  the  model  is  the  stack  voltage  and  the  electrical  power.  The 
thermal  power,  the  stack  temperature  and  the  rejected  heat  are 
computed  under  the  condition  that  the  temperature  of  all  cells  is 
uniform.  The  existing  gas  volumes  in  the  stack  have  a  big  influ¬ 
ence  on  the  performance  especially  for  dynamic  loads,  varying 
gas  supply  or  gas  humidifying.  The  most  important  input  and 
output  parameters  are  listed  in  Table  1 . 

The  gas  diffusion  layer  is  modelled  as  a  porous  media,  which 
is  usually  made  out  of  carbon.  For  the  diffusion  processes  it  is 
necessary  to  calculate  the  amount  of  liquid  water  in  this  media 
because  at  certain  operation  conditions  the  pores  get  filled  with 
water  and  the  diffusion  process  of  the  gases  decreases  (Table  2). 

Also,  the  parameters  and  properties  of  catalyst  layers  and  the 
membrane  are  summarised  in  Tables  3  and  4. 

An  overview  of  all  above-mentioned  different  layers  and  parts 
of  a  fuel  cell  is  shown  in  Fig.  2. 


Table  1 

Stack  parameters 


Stack 

Parameters 

Cell  area 

Number  of  cells  in  the  stack 

Heat  capacity 

Gas  volume  anode 

Gas  volume  cathode 

Initial  conditions  (temperature,  gas  composition,  pressure) 

Input 

Load  current/power 

Cooling  medium  (type,  temperature,  flow  rate) 

Gas  inlet  (pressure,  flow  rate,  temperature,  composition) 

Output 

Voltage/current 

Cooling  medium  (temperature) 

Gas  outlet  (pressure,  flow  rate,  temperature,  composition) 
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Fig.  1.  Fuel  cell  system  model. 


Catalvst  Lavers 

Fig.  2.  Different  layers  within  a  fuel  cell. 


2.2.  Main  relations  incorporated  in  the  dynamic  fuel  cell 
model 


Cell  voltage: 

^cell  =  f(T,  Pi,  A,  jH2,An,  y02,Ca,  yu2,  Ca,  y02,An,  j,  A  Cell,  Cd) 


The  general  dependencies  of  modelled  values  are  given  in 
Section  2.2.1,  then  a  more  detailed  description  of  cell  voltage  is 
given  in  particular  regarding  the  double  layer  capacitance. 


Table  2 


Gas  diffusion  layers 


Gas  diffusion  layer 

Structure 

Porous  media 

Transported  species 

Anode  gases  (H2,  H20_gas,  H20_liq) 

Electrons 

Cathode  gases  (N2,  02,  H20_gas,  H20_liq) 

Phenomena 

Gas  transport  because  of  diffusion  and  convection 
Liquid  and  gaseous  water  transport 

Porosity  depends  on  the  amount  of  liquid  water 

Parameters 

Thickness 

Porosity 

Diffusion  coefficients  for  H20_gas,  H20_liq,  N2, 
02  and  H2 

Cell  temperature: 

T  —  Th?  Tcooi9in,  VcooF  Ccooi,  mCQ n,  stack  size,  ?A ir) 
Table  3 


Catalyst  layers 


Catalyst  layer 

Structure 

Porous  media,  mixture  made  of  the  membrane  and 
backing  layer  material  and  catalyst  particles 

Transported  species 

Anode  gases  (H2,  H20_gas,  H2OJiq) 

Electrons 

Cathode  gases  (N2,  02,  H20_gas,  H20_liq) 

Phenomena 

Double  layer  capacitance 

Electrochemical  reaction 

Water  production  at  the  cathode 

Parameters 

Thickness 

Porosity 

Exchange  current  densities  for  anode  and  cathode 
Diffusion  coefficients  for  H2_gas,  H20_liq,  N2,  02 
and  H2 

Table  4 
Membrane 


Z.  Lemes  et  al.  /  Journal  of  Power  Sources  154  (2006)  386-393 


389 


Membrane 

Structure 

Polymeric,  commonly  made  of  Nation® 

Transported  species 

Protons 

Water 

Permeable  to  hydrogen  and  oxygen 

Phenomena 

Proton  transport 

Water  transport  due  to  diffusion,  hydraulic  pressure 
and  the  drag  effect  of  the  hydrogen  ions 

Swelling 

Proton  conductivity  as  a  function  of  the  membrane 
humidity 

Diffusion  through  the  membrane  H20_liq,  N2,  02 
and  H2 

Parameters 

Thickness 

Number  of  SO3H 

Diffusion  coefficients  for  H20_liq,  N2,  02  and  H2 

Gas  transport  within  the  diffusion  media: 


W,GDL  =  f(T,  Pi,  yir FF,  yi, GDL,  j,  ^ Cell*  G) 


Herein,  £  is  a  function  of  the  amount  of  water  in  the  according 
media. 

Water  transport  within  the  membrane: 


Afif20,Mem  =  f(T,  Pi,  yH20,Ca,  .yH20,An,  j,  A  Cell,  %em) 


Oxygen  transport  within  the  membrane: 


WD2,Mem  —  f(T,  Pi,  yo2,Ca,  J02,An»  j\  ^Cell,  ^Mem) 


Hydrogen  transport  within  the  membrane: 


Afif2,Mem  =  f(T,  Pi,  yH2,Ca,  yH2,An,  J,  A  Cell,  ^Mem) 


Nitrogen  transport  within  the  membrane: 


A^N2,Mem  =  /(T,  Pi,  yN2,Ca,  jN2,An,  j,  ^Mem) 


2.2.7.  Cell  voltage 

The  cell  voltage  is  described  by  the  following  equation: 
t/cell  =  U0-  UAn  -Uq-  f/Ca  (1) 


where  Uq  is  the  open  cell  voltage,  UAn  and  Uca  the  voltage  losses 
on  anode  and  cathode  side  and  Uq  represents  the  ohmic  losses 
due  to  ionic  conductivity  of  the  membrane  and  due  to  electronic 
conductivity  of  bipolar  plates,  diffusion  and  catalyst  layers. 

The  open  cell  voltage  is  a  function  of  temperature  and  partial 
pressures  of  reactants  and  is  calculated  as  follows: 


A  G°  AS 

- + - (T 

IF  2  F 


r°)  + 


RT 

2F 


In 


P  o2 


+  In 


Fig.  3.  Equivalent  circuit  of  a  double  layer. 


Based  on  the  Butler- Volmer-Equation  [6]  which  describes 
the  reduction  and  oxidation  processes  on  the  electrode  and 
neglecting  the  opposite  reaction  on  each  electrode  (Tafel- 
Equation  [6,7]),  the  losses  through  the  charge  transfer  can  be 
described  as  follows: 

RT  j 

Uca  =  - -  ■  In  —  (3) 

<*Ca  •  F  J 0,Ca 

The  exchange  current  density  jo  and  transfer  coefficient  a  are 
measures  for  the  reaction  kinetic.  The  exchange  current  den¬ 
sity  is  concentration  and  temperature  dependent.  Because  of  the 
considerably  higher  exchange  current  density  for  the  hydrogen 
electrode  the  charge  transfer  losses  on  the  anode  side  can  be 
neglected. 

The  ohmic  losses  are  given  by  the  following  equation: 

Uq  =  Re\  •  j  +  7?Mem  •  j  (4) 


where  Rei  represents  the  losses  through  the  electronic  conduc¬ 
tivity  and  Ryiem  is  the  voltage  loss  due  to  ionic  conductivity  of 
the  membrane.  The  conductivity  of  the  membrane  depends  on 
temperature  and  water  content  [8]. 

Regarding  the  fast  load  changes  it  is  of  high  interest  to  know 
the  voltage  response  and  the  storage  capacity  of  a  double  layer. 
Fig.  3  shows  the  equivalent  circuit  of  the  double  layer  and  the 
corresponding  resistances. 

In  this  case  the  double  layer  is  modelled  as  the  parallel 
connection  of  double  layer  capacitance  and  charge  transfer  resis¬ 
tance.  With  Fig.  3  and  Eq.  (3)  following  relationships  are  given: 


RT 

Urd  =  UCd  =  Uca  =  - -  In 

OiCaF 


jcD  =  Gd 


dUcD 

d  t 


Using  Kirchhoff’s  first  law  and  Eqs.  (5)  and  (6)  following  dif¬ 
ferential  equation  for  the  cathode  voltage  is  obtained: 


dGca  j  70, Ca  ( Uca\ 

”d T  ~cj  ' exp  V“^ ) 


with: 

RT 

K  =  - 

UCaF 


Because  of  the  low  charge  transfer  resistance  also  the  result¬ 
ing  time  constant  is  small  so  that  the  double  layer  capacitance 
on  the  anode  side  can  be  neglected. 
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Fig.  4.  Cell  voltage  at  dynamic  load  current  changes. 


Fig.  6.  Simulated  voltage  losses  according  to  dynamic  excitation  in  Fig.  4. 


2.3.  Simulation  results 

In  this  section  the  dynamic  behaviour  of  the  fuel  cell  model 
is  shown  and  compared  with  the  measurements  of  a  stack.  The 
stack  was  air  cooled  and  the  gases  at  the  inlet  were  not  humidified 
during  the  measurement  to  see  the  influence  of  the  membrane 
resistance  during  operation. 

In  Fig.  4  the  cell  voltage  of  the  simulation  is  compared  with 
measured  values  for  dynamic  loads.  There  is  an  excellent  con¬ 
formance  at  dynamic  loads  between  the  simulated  and  measured 
stack  voltage.  The  same  load  characteristics  as  in  Fig.  4  are  based 
on  Figs.  5-7. 

Fig.  5  shows  the  results  of  calculation  of  the  different  over¬ 
potentials  and  the  cell  voltage.  The  cathode  overpotential  has 
similar  characteristic  as  the  current  density  shown  in  Fig.  4. 


The  anode  overpotential  has  similar  characteristic  as  the  cath¬ 
ode  overpotential  but  its  amplitude  is  very  low.  The  membrane 
voltage  loss  depends  on  the  current  density  and  the  proton 
conductivity.  Because  of  the  linear  correlation  between  current 
density  and  membrane  resistance  for  a  current  density  increase 
there  is  a  high  initial  increase  in  the  membrane  overpotential. 
Later  on,  the  resistance  is  decreasing  because  of  the  water  pro¬ 
duction  in  the  stack  and  the  initially  nearly  dry  membrane  is  now 
moistened  (see  also  Fig.  7). 

In  Fig.  6  the  overpotential  for  a  decrease  in  current  density 
from  15  mA  cm-2  to  150  mA  cm-2  is  blown  up.  In  this  case 
the  cathode  overpotential  does  not  increase  by  a  step  because 
of  the  capacitive  double  layer.  The  discharge  of  this  layer  takes 
a  few  seconds  and  during  this  time  the  overpotential  increases 


200  210  220  230  240  250 

Time  /  Min 


Fig.  5.  Simulated  voltage  losses  according  to  dynamic  excitation  in  Fig.  4. 


Fig.  7.  Simulated  membrane  overpotential  and  humidity  according  to  dynamic 
excitation  in  Fig.  4. 
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following  the  characteristics  of  an  exponential  function.  The 
decrease  in  the  oxygen  partial  pressure  results  in  an  overpotential 
increase. 

In  Fig.  7  the  characteristic  amount  of  water  and  the  cor¬ 
responding  membrane  overpotential  is  shown.  The  membrane 
overpotential  is  the  result  of  the  current  density  and  the  proton 
conductivity  which  depends  on  the  water  content.  The  water  con¬ 
tent  is  characterised  as  Lambda,  which  is  defined  as  the  number 
of  water  molecules  divided  by  the  number  of  SO3H  molecules. 
When  the  current  is  low  the  production  of  water  at  the  cath¬ 
ode  is  also  low  and  the  membrane  gets  dry  and  their  resistance 
increases.  Going  to  higher  load  steps  the  overpotential  is  ini¬ 
tially  very  high  because  of  the  dryness  of  the  membrane.  After 
some  time  the  membrane  becomes  more  wet  and  the  resistance 
decreases. 

3.  HIL  test  benches 

For  the  development  of  fuel  cell  stacks  or  fuel  cell  systems, 
there  is  an  increasing  need  in  suitable  test  benches.  The  MAG¬ 
NUM  company  has  been  one  of  the  first  suppliers  of  fuel  cell 
test  benches  and  in  particular  of  HIL  capable  test  benches.  The 
advantage  of  using  HIL  test  stands  is  that  in  early  stage  of  devel¬ 
opment  testing  of  whole  system  can  be  performed,  even  if  all 
the  hardware  is  not  available.  Fig.  8  shows  the  configuration  of 
a  conventional  test  stand. 

As  shown  in  Fig.  8  the  whole  fuel  cell  system  can  only  be 
tested  if  all  sub- systems  and  units  are  available.  In  contrast,  Fig.  9 
shows  the  configuration  if  there  is  a  possibility  to  integrate  HIL 
simulation. 

As  shown  in  Fig.  9  now  it  is  possible  to  replace  defective 
or  not  available  components  with  their  mathematical  models 


Fig.  8.  Conventional  test  bench. 


and  perform  analysis  of  operational  behaviour  of  the  fuel  cell 
system. 

Because  of  the  possibility  to  operate  the  hardware  together 
with  simulation  models  one  of  the  main  aspects  for  the  realisa¬ 
tion  of  HIL  test  stands  is  the  safety.  This  high  safety  standard 
is  reached  through  the  implementation  of  multistage  switch  off 
procedure  on  the  PLC  (Fig.  9).  Further  on,  each  of  the  set  val¬ 
ues  for  the  controllers  calculated  in  the  model  is  checked  for 
the  plausibility  by  the  PLC  before  being  passed  on  to  the  hard¬ 
ware.  The  PLC  will  restrict  the  set  values  to  reasonable  limits 
and  prevent  dangerous  operation  conditions.  This  is  especially 
important  in  case  of,  e.g.  to  mix  hydrogen  and  oxygen  on  the 
anode  side,  or  when  the  pressure  set  values  were  calculated  by 
the  model. 

Another  important  aspect  regarding  automotive  applications 
is  the  dynamic  operation  of  the  test  bench.  Therefore,  fast  mass 
flow  and  pressure  controllers  are  integrated  and  also  fast  data 
acquisition  and  data  exchange  between  the  PLC  and  control  PC, 
respectively,  simulation  PC  are  realised. 

In  order  to  afford  dynamic  changes  of  the  gas  humidification 
in  the  test  stand,  model  based  control  algorithms  are  imple¬ 
mented.  Thermodynamic  simulations  have  been  integrated  to 


Fig.  9.  HIL  capable  test  bench. 
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Fig.  10.  Automation  and  data  transfer  within  HIL  test  bench. 


develop  new  control  algorithms  for  temperature  and  moisture. 
These  methods  allow  to  determine  the  actual  moisture  with  high 
accuracy,  by  measuring  solely  temperature  and  considering  heat 
capacities  of  gases  and  equipment. 

Fig.  10  gives  an  overview  of  the  general  configuration  of  the 
above  described  test  bench  and  the  interaction  of  the  test  stand 
with  the  control  and  simulation  PC. 

As  shown  in  Fig.  1 0  the  fast  data  transfer  between  the  automa¬ 
tion  PC,  HIL  computer  and  the  PLC  is  realised  via  ethernet 
connection.  The  HIL  simulation  can  be  started  by  the  opera¬ 
tor  by  choosing  an  adequate  operation  mode  on  the  automation 
PC.  The  simulation  model  is  running  on  the  HIL  PC  in  the  real 
time  writing  selected  set  values  to  the  PLC  and  simultaneously 
reading  actual  values  from  PLC. 

3.1.  Application  examples 

In  the  following  two  application  examples  of  HIL  simu¬ 
lations  are  presented.  Because  of  the  increasing  interest  of 
automotive  industry  on  the  fuel  cell  technology  there  is  also 
an  increasing  need  for  testing  of  fuel  cell  powered  cars.  An 
example  of  how  the  HIL  capable  test  bench  can  be  used 


Fig.  12.  FuelCell-HIL  simulator. 


in  the  development  process  of  fuel  cell  vehicle  is  shown  in 
Fig.  11. 

In  this  example  the  test  bench  represents  the  power  unit  of 
the  vehicle  where  in  a  fuel  cell  hydrogen  together  with  oxygen 
is  converted  into  electricity,  which  can  be  supplied  to  the  electri¬ 
cal  motor  of  vehicle  propulsion  unit.  On  the  simulation  PC  the 
vehicle  model  is  running  in  real  time.  Using  this  configuration 
(Fig.  11)  static  and  dynamic  behaviour  of  a  power  unit  can  be 
tested  in  an  early  development  stage  when  the  real  vehicle  is  not 
available. 

Next  example  shows  how  the  above  described  stack  model 
running  on  a  real  time  target  (Fig.  12)  can  be  integrated  in  a 
surrounding  system  and  can  be  used  within  HIL  simulation  for 
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Fig.  11.  HIL  simulation  of  a  fuel  cell  vehicle. 
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Fig.  13.  Demonstration  of  a  FuelCell-HIL  simulator. 


design  of  components  and  controls.  In  this  case  the  real  time 
target  system  acts  as  a  stack  dummy. 

The  FuelCell-HIL  simulator  shown  in  Fig.  12  can  be  con¬ 
nected  to  the  surrounding  system  via  the  input  and  output 
ports  and  so  the  system  behaviour  can  be  analysed  and  con¬ 
troller  design  can  be  done  when  the  real  stack  is  not  avail¬ 
able.  An  application  of  the  FuelCell-HIL  simulator  shown  in 
Fig.  13  illustrates  the  interaction  between  the  simulation  and 
the  hardware.  For  the  demonstration  purposes  here  the  load 
current  and  the  pressure  can  be  changed  by  the  user  via  poten¬ 
tiometer  and  the  reaction  of  the  stack  can  be  observed.  On 
the  gauges  the  stack  voltage,  electrical  and  thermal  power  and 
stack  temperature,  which  were  calculated  by  the  model,  are 
shown. 


4.  Summary 

A  dynamic  model  of  PEM  fuel  cell  has  been  presented.  In 
addition  to  the  calculation  of  cell  voltage  this  model  is  also  capa¬ 
ble  of  calculating  time  dependent  temperature  of  the  stack  and 
partial  pressures  of  the  gaseous  species.  The  voltage  losses  on 
the  cathode  are  calculated  considering  the  double  layer  capaci¬ 
tance. 

The  setup  of  HIL  capable  test  benches  has  been  shown  and 
the  special  requirements  for  using  these  test  benches  for  fuel  cell 
testing  have  been  described.  The  use  of  the  HIL  test  benches 
together  with  simulation  models  and  the  resulting  advantages 
for  the  development  of  fuel  cell  systems  has  been  demonstrated. 
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